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ABSTRACT 

We present a detailed analysis of a 3.5 s long bmst from SGR1900+14 which 
occurred on 2001 July 2. The 2-150 kcV time-integrated energy spectrum is well 
described by the sum of two blackbodies whose temperatures are approximately 
4.3 and 9.8 keV. The time-resolved energy spectra are similarly well fit by the sum 
of two blackbodies. The higher temperature blackbody evolves with time in a 
manner consistent with a shrinking emitting surface. The interpretation of these 
results in the context of the magnetar model suggests that the two blackbody 
fit is an approximation of an absorbed, multi-temperature spectrum expected on 
theoretical grounds rather than a physical description of the emission. If this is 
indeed the case, our data provide further evidence for a strong magnetic field, and 
indicate that the entire neutron was radiating during most of the burst duration. 

Subject headings: stars: individual (SGR1900-I-14) - stars: neutron 

1. Introduction 

The soft gamma repeater SGR1900+14 was discovered in 1979 when it emitted 3 short 
bursts of soft gamma-rays in 3 days (Mazets et al. 1979). Its next recorded appearance 
occurred some 13 years later (Kouveliotou et al. 1993). Attempts were made over the 
years to obtain a precise position for the source and identify its counterpart (Hurley et al. 
1994; Vasisht et al. 1994; Hurley et al. 1996), but this remained elusive until 1998, when 
the source entered a new period of activity, allowing it to be localized accurately by the 
interplanetary network (IPN: Hurley et al. 1999a). The precise source location was found to 
be consistent with that of a previously identified ROSAT quiescent X-ray source (Vasisht et 
al. 1994; Hurley et al. 1996). Observations with ASCA further revealed that the source had 
a 5.16 s period (Hurley et al. 1999b), and observations with RXTE demonstrated that the 
period was increasing rapidly (1.1 x 10^^'^ s s^^, Kouveliotou et al. 1999). The counterpart to 
SGR1900-I-14 has not been found yet. If it is associated with the Galactic supernova remnant 
G42.8-I-0.6, it could be as close as 5 kpc (Hurley et al. 1999b). However, the source position 
also appears to be very close to a cluster of high mass stars, and it has been proposed that 
this may be the birthplace of the neutron star (Vrba et al. 2000). If so, its distance could 
be roughly 12-15 kpc. In all this paper we adopt a distance of 10 kpc. 

The 1998 activity of SGR1900+14 culminated in the giant flare of August 27 (Frail et 
al. 1999; Hurley et al. 1999c; Feroci et al. 1999; Mazets et al. 1999), which was followed by 
numerous, smaller bursts (e.g. Ibrahim et al. 2001). The next major period of activity of 
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SGR1900+14 came in 2001 (Guidorzi et al. 2001; Hurley et al. 2001a, b; Feroci et al. 2001; 
Ricker et al. 2001a, b, c; Montanari et al. 2001, Feroci et al. 2004). During this episode, 
it became apparent that this source emits not only the common, short SGR bursts, with 
durations of about 200 ms and fiuences > lO^^crgcm"^, and the much rarer giant flares, 
lasting for minutes and having fiuences > 10^'^crg cm~^, but also, a class of high-fluence 
intermediate bursts, whose durations and fiuences fall somewhere in between, and may in 
fact form a continuum (Kouveliotou et al. 2001; Woods et al. 2003; Feroci et al. 2003). 
Such events had been observed in the aftermath of the 1998 August 27 giant flare, but were 
thought to be related to it. In retrospect, it seems likely that such bursts were also emitted 
by SGR0525-66, in the aftermath of the famous 1979 March 5 giant flare (see Golenetskii et 
al. 1984). The properties of these intermediate bursts are interesting for numerous reasons, 
not the least of which are that some display X-ray afterglows, similar to the one observed 
during the 1998 August 27 event (Feroci et al. 2001; Thompson and Duncan et al. 2001), 
and that their longer duration permits studies of their spectral evolution. 

Duncan and Thompson (1992), Paczyhski (1992), and Thompson and Duncan (1995, 
1996) have proposed that the soft gamma repeaters are magnetars , i.e., neutron stars with 
magnetic fields B ^ lO^^G. In this model, magnetic dissipation causes the neutron star 
crust to fracture, and Alfven waves accelerate electrons, resulting in short (200 ms) bursts of 
soft gamma-radiation. Much more rarely, magnetic reconnection provides the energy for a 
longer, extremely energetic giant flare involving the entire neutron star magnetosphere. The 
periodicity observed in both the quiescent soft X-ray emission from SGR1900-I-14 and in the 
giant flare, the high spin-down rate, and the energetics of the giant flare, are all consistent 
with the main features of the magnetar model. 

In this paper, we analyze an intermediate burst from SGR1900-I-14, which occurred on 
2001 July 2, emphasizing broad band and time-resolved X-ray spectral modeling obtained 
using the data of the FREGATE (FREnch GAmma-ray TElescope) and WXM (Wide Field 
X-ray Monitor) experiments aboard the HETE (High Energy Transient Explorer) spacecraft. 

2. Observations and Instrumentation 
2.1. The FREGATE spectrometer 

2.1.1. Description 

The FREGATE gamma-ray burst (GRB) experiment has been described in detail (see 
Atteia et al. 2003 for the details of the instrument and data modes), so we will only sum- 



-4- 



marize here its most relevant characteristics. It consists of four Nal(Tl) cleaved crystals, 
with total on-axis area 160 cm^, sensitive to photons in the energy range 5-400 keV. It has 
various modes for reading out data, including continuous 128 channel energy spectra with 
5 second resolution, continuous 4 channel energy spectra with 0.16 s resolution, and indi- 
vidually time- and energy-tagged photons which are recorded when a burst trigger occurs. 
This last data type, which allows detailed studies of the spectral evolution of bright bursts, 
consists of 256000 photons (64000 per detector) in 256 energy channels spanning the range 
5-400 keV whose times are tagged to 6.4 //sec. In the present paper, we do not use data 
below 7 keV because of instrumental noise. 

The unique features of FREGATE as far as the present observations are concerned are 
first, its wide energy range, from X-ray to gamma-ray energies, and second, its good energy 
resolution, which ranges from ~ 12% at 122 keV to ~ 42% at 6 keV. The wide energy range 
is particularly interesting for the observations described here, because this is the first time 
that both the X-ray and soft gamma-ray components of an intermediate SGR burst have 
been observed by a single experiment. 

2.1.2. Calibration 

An account of the procedure used to construct the response matrix is given in Olive 
et al. (2003). The spectral response of FREGATE at low energies, which is important for 
the analysis presented here, has been tested in flight with the Crab nebula and with the 
diffuse X-ray background (DXRB). These calibrations led to the discovery of an absorption 
which was larger than expected below 15 keV (e.g. the absorption measured at 8 keV is 32% 
instead of the expected 20 %). This absorption has been attributed to a degradation of the 
surface of the crystal, and corrected empirically in the response matrix. With this correction, 
FREGATE measures the correct parameters for the spectra of the Crab (at angles ranging 
from to 45 degrees) and for the DXRB. In addition, the agreement between FREGATE and 
the WXM for bright GRBs (for which we have enough statistics to do useful comparisons) 
is better than 5% in their common energy range, 7 to 20 keV. We therefore believe that the 
FREGATE spectral response is well known and reliable down to 7 keV at least. 

The spectral analysis of very bright transients may also suffer from dead time and pulse 
pile up effects in the detectors. These effects were carefully measured on the ground before 
launch. The dead time in the burst buffers measured on the ground was found to be ~ 17 fj,s 
(this is the measured dead time resulting from the combination of various electronic dead 
times). For the brightest 25 ms of the burst the dead time was close to 20% , so the true 
counting rate at maximum could have reached 15000 cts/sec per detector above 5 keV. During 
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the decreasing part of the burst, the dead time was between 5% and 10%. The FREGATE 
electronics were specially designed to avoid spectral distorsion at high count rates. Even 
with significant dead time the spectra are not distorted as long as pile up remains negligible. 
Pile up occurs when two photons arrive in a time window of 1.2 fis. This means that pile 
up can affect at most 2% of the counts at the maximum of the peak, at most 1% at the 
beginning of the interval following the peak (interval 14, see below), and at most 0.5% of 
the counts during most of the rest of burst. We thus consider that the spectral analysis 
presented below is not significantly affected by the burst high counting rate. 

2.2. The WXM spectro-imager 

2.2.1. Description 

The WXM is composed of two identical units each consisting of a 1 dimensional coded 
mask and two 1 dimensional position-sensitive proportional counters (PSPCs) placed 187 
mm below the mask. The two units are placed in orthogonal directions. The WXM covers 
the energy range 2-25 kcV, with a maximum effective area reaching 114 cm^ at 8 kcV, when 
the absorption of the mask is taken into account. The WXM has various data readout 
modes, and when a trigger occurs individual photons are recorded, which are tagged in time, 
position and energy. These data were used in section 4 to construct the joint spectrum 
(WXM+FREGATE) of the burst discussed here. 

2.2.2. Calibration 

The energy response of the WXM has been constructed using an analytical model of 
the response measured on the ground. This model takes into account the mask pattern, 
the absorption coefficients for the gas in the PSPCs, the energy-channel conversion and the 
energy resolution measured within the volume of the detector. The response of the WXM 
has been tested in flight with the Crab nebula, measured at angles varying from 0° to 30° 
off-axis. The spectral slope is correctly measured, as well as nn, the column density of the 
gas, except at large incident angles 30°) where uh appears to be underestimated. This 
last feature may be relevant for the burst studied here, which arrived at an angle of 32.7°. For 
bursts with count rates above ~ 4000 c/s, the WXM experiences significant dead time. In 
the present case. The dead time in the WXM varies significantly during the burst, resulting 
in a WXM light curve which does not closely follow the light curve recorded by FREGATE 
at low energies (see Fig. 2 of Torii et al. 2003). In addition the housekeeping data needed 
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to measure the dead time of the WXM were only available on a 5 s timescale for this burst. 
These two factors have the consequence that we can use the WXM data only for the study 
of the time averaged spectrum of the July 2 burst (Section 4). The time- resolved analysis 
(Section 6) must be based on FREGATE observations only. 

2.3. Context: the 2001 activation 

The 2001 activation of SGR1900-I-14 began on April 18 and ended around July 8. During 
that period it emitted over 100 bursts. The exact number is uncertain for two reasons. First, 
many events were detected by a single instrument in the IPN and could not be localized. 
Second, the source is known to emit weak and/or soft spectrum events (Woods ct al. 2001, 
2003), and these would have been undetectable to the experiments in the IPN. However, 
among the events whose origin is definitely SGR1900+14, two are noteworthy for their long 
durations and high fluences: that of April 18 (Guidorzi et al. 2001, Hurley et al., 2001a, 
Feroci et al. 2003), and the present 2001 July 2 event. 

This burst (HETE trigger 1576) was detected by the FREGATE experiment at 12846.529 
s UT and also by the Wide Field X-Ray Monitor aboard HETE, and both automated and 
ground-processed locahzations were sent out. The WXM error circle, which has a radius of 
12 ', includes the position of SGR1900-I-14 (Ricker et al. 2001c). In addition, this event was 
observed by a number of spacecraft in the IPN: Ulysses, Mars Odyssey (HEND experiment) 
BeppoSAX (GRBM experiment), and Wind (Konus experiment). The triangulated position 
is also in agreement with that of SGR1900+14(Hurley et al. 2001a). At the time of the 
burst, the HETE detector axis was pointed at «(2000) = 280.4°, 5(2000) = -22.8°, 32.7 
degrees from SGR1900+14. Thus the effective geometrical area exposed to SGR 1900+14 
was 135 cm^ for the detectors of FREGATE and 14.6 cm^ for the WXM detectors. 



3. Temporal analysis 

The 7-100 keV time history of the 2001 July 2 burst from the time-tagged FREGATE 
data is shown in figure 1 where the reference time is the HETE trigger time {Tfr). FREGATE 
recorded a total of 1.07 x 10^ counts from the burst. At maximum the total count rate on a 
25 ms timescale reached 43000 cts/sec. 

Using the continuous 0.16 s, 4 channel FREGATE data and the WXM data, we have 
searched for evidence of SGR activity before (i.e. precursors) and after (i.e. extended tail) 
the main peak which triggered FREGATE. With the possible exception of a very weak tail. 
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discussed briefly in section 4, no such activity was found. The sensitivity of FREGATE for 
a 1 second long signal with an SGR spectrum is of the order of 4 x 10~^ erg cm~^ s~^. This 
sensitivity guarantees the detection of a precursor 1.5 times fainter than the one observed by 
Ibrahim et al. (2001) before the burst of August 29th. It is however, much too low to permit 
the detection of a tail which would be comparable to the tail of the August 29th burst. The 
instrumental background, plotted as a dotted line in figure 1, was interpolated using the 
count rates prior to T^^ — 2 s and after Tfr + 6 s. The total duration of this intermediate 
burst is Tgo = 3470 ms; this may be compared to a typical SGR burst duration of ~ 200 ms. 

Motivated by the light curve morphology, wc divided the burst into 9 consecutive and 
contiguous intervals (hereafter designated II to 19). The burst begins with a sudden rise in 
the count rate (II, lasting 45 ms, with an e-folding rise time of 32 ms), followed by an initial 
spike (12, with a duration of 30 ms) and a plateau (13, with a duration of 70 ms). This is 
followed by a smooth decay with an e-folding time of 3.9 s, lasting ~ 3400 ms, that we have 
divided into 4 equally long intervals (14 to 17), and a final, faster decay lasting 300 ms (18, 
duration 380 ms)^. The event ends with a faint ~ 2s long decreasing tail (19). 



4. The 2-150 keV time-integrated spectrum 

4.1. Spectral modeling 

We have used XSPEC (Arnaud 1996) to deconvolve the WXM and FREGATE count 
spectra. In all our spectral analyses wc have taken into account a Galactic absorption tlh = 
2.5 X 10^^ cm^^, determined from the ASCA measurements of the quiescent X-ray counterpart 
(Hurley et al. 1999b). Since there are no source counts in the spectrum above 150 keV, we 
first attempted to fit the 2-150 keV time-integrated count spectrum (taken over II to 18) 
using several single component models. We tried a single temperature optically thin thermal 
bremsstrahlung spectrum (OTTB), a simple or broken power law (PL, BKPL) and a single 
blackbody (BB). None of these models reproduces the shape of the measured spectrum from 
2 to 150 keV, and all of them may be rejected on statistical grounds (x^/dof= 2146/95, 
19135/95, 890/93, 1965/95 respectively). A similar conclusion applies to the time resolved 
spectra of intervals II to 19 taken separately. The OTTB function has been successfully 
used to fit the spectra of many of the short SGR bursts above 15 keV. In the present case 
it fails to reproduce the rollover of the spectrum observed by FREGATE below 15 keV (an 



^Fig. 2 of Olive et al. (2003) was erroneous. It showed a spike just before the decay which was an artefact 
due to an uncorrected timing error in the time-tagged photon data. 
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illustration of this failure can be found in Fig. 3 of Olive et al. 2003). We note that a similar 

rollover in the spectrum of bright SGR bursts, was previously reported by Fenimore et al. 
(1994), based on observations of SGR 1806-20 with the International Cometary Explorer. 
More recently, Feroci et al. (2004) also reported a deficit of photons below 15 keV in 10 
short bursts from SGR 1900+14 observed with the BeppoSAX narrow field instruments. 

As an alternative, we have tried the magnetar spectrum of Perna et al. (2001). In 
this model, the spectrum originates as thermal emission from the neutron star surface, but 
includes the effects of anisotropic heat fiow through a magnetized neutron star envelope, 
reprocessing in a light-clement atmosphere, and rclativistic corrections. R. Perna has kindly 
extended her results to the higher energy range of FREGATE to make this possible (private 
communication, 2002). However, we do not obtain an improved fit using it. 

Next we tried several two-component models, which included combinations of PL, 
BKPL, and OTTB. They also failed to reproduce the shape of the spectrum: the best 
fit was given by a BKPL-I-BKPL model, which is however clearly rejected (x^/dof= 358/89). 
We eventually found that the only two-component model that can be fit to the experimen- 
tal spectrum is composed of the sum of two blackbodies (BB) with different temperatures. 
Hereafter we will refer to these two thermal components as the Low and the High Tem- 
perature Blackbodies (LTBB and HTBB respectively). The experimental spectrum and its 
residuals are plotted in figure 2. The temperatures for this time-integrated spectrum are 
kTLTBB = 4.25 ± 0.1 keV and kTuTBB = 9-8 ± 0.3 keV. 

The reduced value that we derive for this fit (x^/dof= 137/93, corresponding to a null 

hypothesis probability of 0.002) is still not completely satisfactory. However, the value can 
be reduced to x^/dof= 104/93 (null hypothesis probability = 0.2) by adding a 3% systematic 
error in quadrature to the statistical errors to account for calibration uncertainties affecting 
this high signal spectrum. Thus we conclude that the HTBB+LTBB model provides an 
acceptable fit to the measured WXM+FREGATE spectrum over the 2 to 150 keV energy 
range. 

4.2. Two vs. three blackbodies 

Since the temperatures of the two blackbodies in our fit differ by not much more than a 
factor of two, the two blackbodies are not detected individually, and this raises the possibility 
that the data could favor a model based on a sum of blackbodies with a narrow distribution 
of temperatures. To answer this question we tried to fit the time-integrated spectrum with 
a three-blackbody model (the three BB model intending to simulate a narrow distribution 
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of temperatures) . 

The three-blackbody model gives a good fit with temperatures equal to 1.4, 4.4, and 10.0 
keV, and corresponding radii of 55, 23, and 4.0 km (at 10 kpc). The fit is fully acceptable 
(x^ = 81 for 91 dof), and provides some improvement over the two-blackbody model (x^ = 
104 for 93 dof). A close examination of the three BB fit shows that the additional blackbody 

component (with kT=1.4 keV) is used to fit the few points below 5 keV which are above the 
two BB fit (see figure 2). Since the spectral response of the WXM at large incident angles 
is subject to larger uncertainties (section 2.2.2) and since the two BB fit is statistically 
acceptable (section 4), wc conclude that, to within the limits imposed by the data, a third 
blackbody component is not strongly required, and we restrict the following discussion to 
the two-blackbody model. 

4.3. Burst energetics 

The spectral analysis was used to derive the burst energetics. A 2-150 keV fiuence of 
19 X 10~^ erg cm~^ and a ~ 3.5 s duration place this event firmly in the "intermediate" 
class^. Specifically, this event looks similar to other intermediate bursts from SGR 1900-1-14 
described in the literature, such as 981028b in Aptekar et al. (2001) (duration = 4 s, fluence= 
48 X 10~^ erg cm~^), and the "unusual burst" extensively discussed by Ibrahim et al. (2001) 
which had a duration of 3.5 seconds and a fiuence above 25 keV of 19 x 10~^ erg cm~^ (vs. 
7.8 X 10~^ erg cm~^ for the event discussed here). 

Assuming for the sake of concreteness a distance of 10 kpc, the average luminosity in 
the 2-150 keV energy range is ~ 6.0 x lO^'^ erg s~^ and the energy release is 2.1 x 10^^ erg. 

5. Compcirison with previous work 

5.1. Spectral analyses of intermediate bursts 

Detailed spectral analyses are available for two other intermediate bursts from SGR1900-I-14: 
the 1998 August 29 event observed with BATSE and RXTE, and the 2001 April 18 event 
observed with the BeppoSAX Gamma- Ray Burst Monitor (GRBM). It should be noted that 
FREGATE covers energies below 25 keV, which were not covered by BATSE or the GRBM, 



^Thc values given here arc the ones measured with FREGATE. Note that, considering the amount of 
dead time (see section 2.1.2), the fiuence given here should be increased by 10% to obtain the true fiuence. 
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but that the FREGATE data do not constrain the burst spectrum above 300 keV, due to 
decreasing sensitivity at these energies. 

Ibrahim et al. (2001) found that the spectrum of the 2001 August 29 burst was correctly 
fitted with an OTTB (kT=20.6 keV) above 25 keV. We find that an OTTB is also a reasonable 
fit to the 2001 July 2 burst spectrum above 25 keV (xVdof = 57.2/47) with kT = 23.0l^j^ 

kcV. Thus the results reported in this paper are consistent with those of Ibrahim et al. 
(2001); indeed, it is possible that, had <25 keV data for that event been available, a two 
blackbody model might have been preferred for it, too. 

Overall, the burst of 2001 April 18 is longer (~ 40 s) and more energetic (with a 40-300 
keV fluence ~ 120 x 10~^ erg cm~^) than the event discussed in this paper. Guidorzi et 
al. find that the spectrum of their bright burst is well fitted by a combination of three 
spectral components (see their Table 2): a 14 keV blackbody, a broken power law with 
a break at 73 keV, and a single power law with an index close to —0.5. Replacing the 
blackbody by a 32 keV OTTB gives an equally good fit. The single power law dominates 
the spectrum at energies above 300 keV, but we will not discuss it here because it is outside 
the FREGATE energy range. We tried to fit the spectrum of our event with a model which 
included a blackbody and a broken power law. The best fit is obtained with a blackbody of 
temperature 5.1 ±0.2 keV plus a broken power law with a break at 53.5 ±5 keV. However, 
the reduced is 1-75 (x^= 112 for 64 dof), indicating that this fit is not acceptable. While 
the FREGATE burst cannot be fit with the sum of a blackbody and a broken power law, 
Guidorzi et al. mention that a fit with the sum of two blackbodies cannot be ruled out for 
the burst of 2001 April 18. C. Guidorzi kindly provided us the best fit parameters of their 
two blackbody fit in the energy range 40-300 keV (private communication, 2003). Their 
temperatures are 13.6 (-0.2;+0.1) keV and 35 (-3;+2) keV, with 72.6 for 66 dof. While 
their low temperature blackbody could be compared with our high temperature blackbody, 
we found no evidence in the FREGATE data for a 35 keV blackbody. We tried to fit our 
40-300 keV spectrum (corresponding to the energy range of the Guidorzi et al. analysis) 
with two EEs with temperatures of 13 keV and 35 keV, but found that convergence was 
attained for two blackbodies having temperatures of 4.1 keV and 10.6 keV (x^= 53.0 for 66 
dof). If we freeze the temperatures of the EEs to the values obtained by Guidorzi et al. we 
obtain a of 568 (for 68 dof), and the 35 keV BB has a normahzation of 10~^, indicating 
that it is rejected by the data. We therefore conclude that the intermediate bursts of 2001 
July 2 and 2001 April 18 have very different energy spectra. 
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5.2. Spectral analyses of short bursts 

In parallel with the analysis of the intermediate burst reported here, we have performed 
the spectral analysis of several classical short bursts. We find that the two BB model 
adequately fits the spectrum of short SGR bursts with parameters similar to those given 
in Table 1. A more detailed discussion of the spectral modeling of short SGR bursts is in 
progress and will be reported in a future paper ( Atteia et al. in preparation) . Recently Feroci 
et al. (2004) presented broad band spectra of ten short bursts from SGR 1900+14 observed 
with the narrow field instruments of BeppoSAX. They find that the OTTB model, which 
successfully fits the spectrum of short SGR bursts above 15 keV, is not acceptable below 15 
keV. They also find that a two BB or a cutoff powerlaw model provide good spectral fits to 
the short SGR bursts in the energy range 1.5 to 100 keV. With the two BB model, the ten 
bursts are jointly fit with a model having temperatures at 3.3 keV and 9.5 keV, similar to 
the temperatures reported in our Table 1. 

Following Feroci et al. (2004), we tried to fit the WXM+FREGATE spectrum with a 
cutoff powerlaw model and find it to be excluded by the data (x^ = 170 for 94 dof, with 
3% systematic error, vs. 104 for 93 dof for the two BB model). The corresponding null 
hypothesis probability is 2 x 10~^ vs. 0.2 for the two BB model. 

6. Time- resolved analysis 

In the first part of this section we display the spectra of intervals II to 19, illustrating the 
validity of the two BB model on short timescales (section 6.1). In a second part (section 6.2) 
we study the evolution of the model parameters in 24 sub-intervals, chosen to provide the 
finest temporal resolution while keeping enough counts to perform a statistically significant 
analysis. 

6.1. Time-resolved spectra 

To approximate the instantaneous spectral shape, we have fit the spectra of intervals 
II to 19 (figure 1) with a two BB model. As explained in section 2.2.2, this study can only 
be done with FREGATE, and therefore it is restricted to the energy range 7-150 keV. We 
successfully fit the spectra of II to 18 with the two BB model, while 19 could be fit with a 
single BB component. Table 1 summarizes the results. Here kTLXBB and kTnTBB are the 
temperatures of the two blackbodies, Rltbb and Rhtbb are the radii of the emitting regions, 
and Lltbb and Lhtbb are bolometric luminosities in units of 10^° erg s~^ for an assumed 
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distance of 10 kpc. In figure 3, we have plotted the 9 unfolded spectra and compared them 
to the model (the spectrum of 19 was fit with a single blackbody model). For all spectra, 
the fit is fully acceptable, the worst case being for 13 (x^/dof = 81.5/66, corresponding to a 
null hypothesis probability of 0.1). In figure 4, we have plotted the 9 count spectra and the 
residuals of the fit. Table 1 and figure 3 show that i) the two blackbody fit works on short 
timescales, and ii) the temperatures of the two blackbodies and the overall spectral shape 
evolve significantly during the burst. 



6.2. Temporal evolution of the spectral peirameters 

In order to study accurately the variations of the spectral parameters along the burst, 
we have sub-divided this event into 24 intervals, noted il to 124 and corresponding to the 
former II, 12, 13, 14 to 17 divided into 19 sub-intervals, 18, and 19. These intervals are the 
shortest ones that can be selected with good statistics. 

We have applied the fitting procedure with the LTBB+HTBB model to the resulting 
24 spectra. The derived parameters (blackbody radius and temperature) arc plotted as a 
function of time in figure 5 and reported in Table 2. In figure 5, the dashed horizontal 
lines indicate the value of the parameters measured for the time-integrated spectrum. The 
luminosities of the two components as a function of time are plotted in figure 6. A crude 
analysis of the variability of the four spectral parameters shows that the two temperatures, 
as well as the radius of the low temperature blackbody, display comparable fiuctuations : 
(t(/o^(Tltbb)) = 0.068, ^(^(Thtbb)) = 0.061, and cT(%(i?LTBB)) = 0.068. There is no 
statistically significant evidence that either temperature or Rltbb evolves with time. The 
radius of the high temperature blackbody, on the other hand, shows significantly larger 
fiuctuations, with cr(/o5'(i?HTBB)) = 0.197, and it tends to decrease with time. 

In order to test for possible correlations between the fit parameters we have plotted 
in figure 7 the temperature as a function of the radius for the two blackbodies, for the 24 
intervals except intervals il, 123, and 124 whose parameters are not well constrained. The 
two blackbodies exhibit different emission patterns: while the radius of the low temperature 
blackbody is nearly constant and independent of the temperature, the high temperature 
blackbody shows a clear anti-correlation between the radius and the temperature. A cor- 
relation analysis using the Pearson test shows that this anticorrelation is significant at the 
4.4 a level. The shapes of the ellipse errors, however, indicate that the observed correlation 
is at least partly due to the fitting procedure. Indeed, the parameter relationship in the 
blackbody model is strongly non-linear. 
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7. Discussion 

7.1. Temperatures, radii, and luminosities 

A first look at the parameters given in Table 1 calls for the following remarks. The 
temperature of the HTBB is very close to the typical temperature of a fireball trapped at 
the surface of a magnetar (~11 keV) predicted hy Thompson & Duncan (2001). The radii of 
the HTBB and LTBB are consistent with those of neutron stars and/or their magnetospheres. 
For distances less than or greater than the assumed 10 kpc, the radii would scale linearly 
with the distance. Finally, the time-integrated luminosities \jltbb and Lhtbb emitted by 
the low temperature and high temperature components differ by less than 20%.^ 

7.2. Comments on the spectral modeling 

We have shown that a two BB model fits the data. This model can be interpreted 
either as a physical or a phenomenological description of the data. In the first case, the 
model, taken literally, implies that we are observing two distinct radiating volumes, and the 
measured parameters give us insight into the physics of these regions. In the second case, 
the fit is only a convenient representation of the emission and little or no physical meaning 
can be attributed to the measured parameters. At this stage, we believe that the study of a 
single event is not sufficient to reach a firm conclusion about the physical meaning of the 2 
BB description. We expect that the study of a large number of SGR bursts will eventually 
resolve this uncertainty, and in the meantime we discuss below the implications of these two 
possible interpretations. 

7.3. The two blackbody interpretation 

7.3.1. The HTBB: A trapped fireball ? 

Two observations lead us to explore here the hypothesis that the HTBB originates from 
a fireball trapped on the surface of the neutron star, as described in Thompson & Duncan 
(2001). The first is the remarkable agreement of the temperature of the HTBB inferred from 
the FREGATE data with the temperature predicted theoretically for a trapped fireball (9.8 



■'While the integrated luminosities are comparable, figure 7 and table 1 show that the quasi-equality of 
the luminosities of the two blackbodies does not hold on short timescales, particularly at early times. 
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vs. 11 keV). The second is the contracting radius of the HTBB blackbody which is invoked 
in the trapped fireball model. 

Can a trapped fireball explain the temporal evolution of the HTBB emission? The 
abrupt disappearance of the emission in about 0.3 s (figure 1, interval 18) is too rapid to be 
due to the rotation of the neutron star, but it is consistent with the behavior of a trapped 
fireball, and we believe that it reflects the intrinsic evolution of the emission. To check 
this hypothesis we tried to fit the light curve above 40 keV, where the HTBB is dominant, 
with a trapped fireball model following the temporal evolution described in Thompson & 
Duncan (2001). The light curve and the fit arc displayed in figure 8. The light curve exhibits 
two episodes of fast decline which cannot be fit with a single trapped fireball. It is nicely 
fit, however, with the sum of two trapped fireballs having evaporation times of 1.3 sec and 
3.8 seconds, and dimensions x — 0-4 and x = 0.1 (see Thompson & Duncan 2001). For 
comparison, x = 2 for a spherical fireball, and x ~ 4 for the 1998 August 27 event (Feroci 
et al. 2001). 

We also note that we do not observe any modulation of the emission with the 5.16 s 
period of SGR1900+14, even though the burst is significantly longer than half the rota- 
tion period. The fact that the 5.16 sec modulation has been observed many times (in the 
quiescent soft X-ray emission and during the giant fiarc), indicates that the rotation axis 
of the magnetar is not directed towards the Earth. It is still possible, however, to have it 
not too far from the line of sight, allowing large regions to remain permanently visible and 
unmodulated. If this burst took place in the polar region continuously visible from Earth, 
this would explain the lack of rotational modulation. 

7.3.2. The LTBB: thermal emission from the heated neutron star 9 

While the energy radiated by the two blackbodics is roughly equal, their radii and 
temperatures are clearly different. The size of the region emitting the LTBB is about 30 
times larger than the average size of the "fireball". One possible explanation is that part 
of the fireball energy is reprocessed in a larger region and emitted at a lower temperature. 
This would also explain the quick drop in luminosity of the LTBB, a characteristic which 
it shares with the HTBB. Under this hypothesis, the stability of the radius of the LTBB 
blackbody shows that the "reprocessing region" is not strongly affected by the amount of 
energy coming from the fireball. 

Could the LTBB be emitted by the surface of the neutron star? We check below whether 
the size of the region emitting the LTBB is compatible with the typical surface area of a 
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neutron star. Under this hypothesis, the LTBB spectral parameters (equivalent radius and 
temperature) are affected by the gravitational redshift of the neutron star. For a given 
source distance, one can derive a relation between the true radius and the gravitational mass 
of the star. These relations are plotted in figure 10. For example, at a distance of 10 kpc, 
a neutron star with a mass of 1.4 M© has a radius of 23 km. Distances of 4, 6, 8, 12, and 
14 kpc give radii of 5, 12.5, 18, 29, and 34 km respectively. Given the uncertainty in the 
distance of SGR 1900+14 (see Section 1), our observations do not exclude the possibihty 
that the LTBB is emitted by the surface of the neutron star, but the numbers above clearly 
show that the whole surface must radiate, a fact which could help explain the absence of 
rotational modulation. 



7.4. The 'magnetcir spectrum' interpretation 

Radiation transfer in the strong magnetic fields of magnetars has been studied in some 
detail in recent years. Here we discuss the interpretation of our data in terms of those 
theoretical models of SGR bursts which take into account the detailed radiation transfer in 
superstrong magnetic fields. This is mainly based on the work of Ulmer (1994) and Lyubarsky 
(2002) which describes the expected X-ray spectrum of short SGR bursts. These authors 
predict that the photosphere of the fireball should radiate a blackbody-like spectrum with a 
typical temperature of ~ 10 keV. However, they note that the opacity of a fireball in a strong 
magnetic field is a function of energy, allowing the observer to see deeper into the fireball at 
lower energies. This has the double efi^ect of broadening the spectral distribution (extending 
it towards lower energies), and lowering the peak of the spectrum. Magnetic photon splitting 
in the strong magnetic field of magnetars further complicates the observed spectrum (Baring 
1995, Lyubarsky 2002). This effect, like double Compton scattering, provides a source of 
additional photons, allowing a very efficient thermalization of the fireball. It also modifies 
the shape of the radiated spectrum. Given the importance of this process, it is interesting to 
investigate whether it can degrade the energy of the photons from a trapped fireball down to 
a few kcV. This issue has been studied by Baring (1995), and also by Thompson & Duncan 
(1995), who find that it can degrade the energy of the photons down to 10 keV, but probably 
not below. This is however a very complex issue, involving both the details of the radiation 
processes and the geometry of the emitting region, and it deserves additional investigation 
in view of the present observations. 

Although the present event is not a short burst, we checked the agreement of our data 
with these theoretical predictions by comparing the WXM+FREGATE spectrum with the 
spectrum predicted by Lyubarsky. Figure 9 shows a good agreement between the observa- 
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tions and the theory below 30 keV. Above 30 keV, the model under-predicts the observed 
counts, even though the effects of photon splitting have been included. The addition of a 
10 keV blackbody as a second component would rectify this (not shown). We can use the 
spectral fit below 30 kcV to infer the effective temperature Tcff of the fireball. We find 
an effective temperature of 6 kcV. This value allows us to compute the size of the emitting 
region, which is A = L/aT^^^, assuming isotropic emission. Taking L = 6 10^° erg cm~^ 
s~^, and Te// = 6 keV, gives A = 4500 km^, which corresponds to the surface of a sphere of 
radius 18.9 km, assuming a distance of 10 kpc, and ignoring the effect of the gravitational 
redshift. 

The fact that this value is close to the radii of known neutron stars further reinforces 
this interpretation. Figure 10 shows the mass-radius relation of the magnetar in this inter- 
pretation. We finally note that the observed luminosity of the burst (L = 6 10^*^ erg cm~^ 
s~^) agrees well with the magnetic Eddington luminosity discussed by Ulmer (1994) and by 
Thompson & Duncan (1995, Section 3.1). Despite the good agreement with theoretical pre- 
dictions, we believe that it is essential to study more bursts, emitted by several sources, to 
definitely validate this interpretation of the data. 

8. Conclusions 

We have presented HETE-FREGATE and WXM observations of an intermediate burst 
from SGR1900-I-14. This is the first time that time-resolved spectral analysis of such a burst 
has been possible over a wide energy range, with good spectral and temporal resolution. We 
have found that the spectrum is well described by the sum of two blackbodies. One possible 
interpretation is that we are observing two distinct emitting volumes. Another is that we are 
observing emission which is affected by radiation transfer effects in a superstrong magnetic 
field. Both interpretations are consistent with the basic features of the magnetar model of 
soft gamma repeaters. However, it is not possible to choose between the two on the basis 
of a single observation. Fortunately, HETE-2 has observed numerous SGR bursts from both 
SGR1900-I-14 and SGR1806-20. Analysis of these observations is underway, and may provide 
less ambiguous support for one of these hypotheses. 

We are grateful to R. Perna for extending her magnetar spectrum model to the FRE- 
GATE energy range and making the results available to us, to C. Guidorzi for giving us 
detailed information about the BeppoSAX data, and to M. Baring for useful comments. KH 
is grateful for HETE support under MIT contract MIT-SC-R-293291. 
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Fig. 1.— The 7-100 keV time history of the 2001 July 2 burst from SGR1900+14 as observed 
by FREGATE (time resolution 20 ms). Inset: initial part of the burst (time resolution 5 
ms). The interpolated background is indicated by a dashed line. The start time corresponds 
to 12846.529 UTC (HETE trigger time). Horizontal lines have been drawn and labeled to 
indicate the intervals corresponding to the 9 spectra of figures 3 and 4 (II to 19). For the 
detailed time-resolved spectral analysis, intervals 14 to 17 have been further subdivided into 
19 sub-intervals (not shown). 
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Fig. 2. — Top: Deconvolved (photon) spectrum of the time-integrated burst (II to 18, see 
figure 1) using a two blackbody model. The two blackbodies and their their sum are plotted. 
The filled circles correspond to the FREGATE 7-150 keV spectrum and the open circles to 
the WXM 2-25 keV spectrum. Middle: Observed (count) spectrum of the time-integrated 
burst and the best fit model (dotted line). Bottom: Residuals between the data and the 
model normalized to the standard deviation. 
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Fig. 3. — FREGATE deconvolved (photon) spectra and the best fit model (the sum of 2 
blackbodies, shown as dashed hnes, and their sum, shown as a sohd hue) for the 9 intervals 
of the burst of 2001 July 2 indicated in Figure 1 (II to 19). Interval 19 can be fit with a 
single blackbody. The spectral parameters for these 9 spectra are given in Table 1. 
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Fig. 4. — FREGATE count spectra and residuals for the two blackbody model for the 9 
intervals of the burst of 2001 July 2 indicated in figure 1. The two horizontal dashed fines 
indicate ± 1 sigma confidence limits on the residuals. The last interval can be fit with a 
single blackbody. The spectral parameters are given in Table 1. 
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Fig. 5. — Spectral parameters of the two BB best fit as a function of time for 24 time intervals 
il to 124 (see text). The bottom panel displays the temperatures of the LTBB (filled circles) 
and HTBB (open circles). The middle panel displays the radii of the LTBB (filled circles) 
and HTBB (open circles), for a source distance of 10 kpc. Note the decreasing trend of the 
radius of the HTBB The error bars arc for 90% confidence level (defined as the parameter 
interval for which the is less than Xmin + 2-7 for variation of a single parameter, re- fitting 
the spectrum at each step). The upper panel shows the null hypothesis probability of the 
fit, and indicates that all the fits are acceptable. 
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Fig. 6. — Bolometric luminosities of tiic LTBB and HTBB as a function of time, for tiie 24 
fits (il to i24, see Table 2) 
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Fig. 7. — Left panel: The radii of the two blackbodies as a function of their temperatures, 
for intervals 12 to 122 (see Table 2). Filled symbols are used for the LTBB and open symbols 
for the HTBB. The circles show the parameters for intervals 12 to ill, and the crosses for 
intervals 112 to 122. Right panel: Contour plot of the 90 % confidence region ellipses for 
each data point. Note that the radius of the LTBB (filled symbols) does not vary with 
the temperature, while the radius of the HTBB (open symbols) is anti-correlated with the 
temperature. 
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Fig. 8. — Fit of the light curve above 40 keV with the sum of two trapped fireballs having 
evaporating times of 1.3 and 3.8 seconds and dimensions x = 0.4 and x — O-l- The time 
resolution is 20 ms. The first 0.25 seconds of the burst (the first spike) are not included in 
the fit. 
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Fig. 9. — Comparison between the modified blackbody spectrum proposed by Lyubarsky 
(2002), with an efi^ective temperature of 6 keV, and the observed spectrum. Although photon 
splitting has been taken into account, the observed spectrum only fits the modified spectrum 
below 30 keV. 
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Fig. 10. — Magnetar mass-radius relation for various distances and models. The dotted line 
is for the two blackbody (LTBB + HTBB) hypothesis. The solid line is for the modified 
blackbody spectrum of Lybarsky (2002) with an effective temperature of 6 keV. The hori- 
zontal dashed hue indicates a mass of 1.4 Mq. There are no solutions in the region to the left 
of the other dashed line. Solutions exist for source distances of ~ 6 — 10 kpc with inferred 
radii which are compatible with the accepted values of neutron star radii (~ 12 — 18 km) 



-30- 



Table 1. Parameters of the two-blackbody fits. Intervals II to 19 correspond to those of 
Figure 1. The columns contain the name of the interval, the duration of the interval, the 
LTBB temperature, radius and bolometric luminosity, the HTBB temperature, radius and 

bolomctric luminosity, the number of degrees of freedom (dof) of the fit, and the Chi 
square per dof. The subscript LTBB (HTBB) corresponds to the Low (High) Temperature 
Blackbody component. The error bars are for the 90% confidence level. 



Interval 


Duration 






^LTBB 


kT jfTBS 


^HTBB 




dof 


xVdof 




sec 


keV 


km 


10^° erg 


keV 


km 


10^° erg 






11-18'^ 


4.093 


4.25t°;i 




2.92 ±0.09 




4.69t°1 


2.59 ±0.09 


99 


1.12= 


11-18'' 


4.093 


4.31^0;} 




3.00 ±0.09 




4.54+0;^ 


2.56 ±0.09 


66 


0.83^= 


lib 


0.045 






2.22 ±0.56 




10.39+43 


5.01 ±0.71 


66 


1.03 


12^ 


0.030 




28.17^21 


14.38 ±1.59 




2.61+t? 


3.42 ±1.53 


66 


0.50 




0.073 


5.49t°;l 


29.17i|2 


10.08 ±2.13 


11.07+ig 


4.58t^° 


4.09 ±2.17 


66 


1.23 


14b 


0.891 


4.45t°;^ 




4.13 ±0.29 


9.27t°1 


7.07t°;^ 


4.79 ±0.31 


66 


1.15 


15b 


0.892 


4.29toj 


26.89ti:3 


3.18 ±0.17 


9.70toj 


4.42t^-j 


2.25 ±0.18 


66 


0.94 


16^ 


0.891 


4.05t°;i 


27.02ti:^ 


2.56 ±0.12 


9.87t°oi 


q oo+0'7 
o.88_o.6 


1.86 ±0.13 


66 


0.99 




0.892 


3.86^0-2 


26.36+?? 


2.01 ±0.09 


10.17^E5;6 


3.52t°;i 


1.72 ±0.10 


66 


0.82 


18^ 


0.379 


3.68l°J 




0.76 ±0.21 


7.71+\i 


3.92^2-^ 


0.70 ± 0.23 


66 


1.22 


19b 


1.94 






6.23t°i 


3 1 o+o> 


0.19 ±0.03 


68 


0.72 



Joint spectrum (WXM and FREGATE) in the energy range 2-150 keV. See Figure 2 
^FREGATE-only spectrum in the energy range 7-150 keV. 
'^This spectrum includes 3% systematic errors. 
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Table 2. Parameters of the two-blackbody fits for the 24 sub-intervals il to 124 (sec text). 
The columns contain the name of the interval, the duration of the interval, the LTBB 
temperature, radius and bolometric luminosity, the HTBB temperature, radius and 
bolometric luminosity, the number of degrees of freedom (dof) of the fit, and the Chi 
square per dof. The subscript LTBB (HTBB) corresponds to the Low (High) Temperature 
BlackBody component. The error bars are for the 90% confidence level. 
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32.42t^;^ 


2.61 ± 0.43 


8.8ltJ-i 


7 "^1 +21 
' •"^-'--2.0 


4.19 ±0.46 


66 


1.18 


il2 


0.158 


4.86t°i 


24.93li;g 


4.49 ±0.56 


11.81+2-4 


2.63+1-^ 


1.75 ±0, 


.56 


66 


0.97 


il3 


0.184 


4.0710-3 


29.02l|^ 


3.01 ±0.57 


8.84tJ;0 


5.5lt};^5 


2.41 ±0, 


.59 


66 


0.80 


il4 


0.199 




27.3i;^:g 


2.79 ±0.62 


8.39il;^ 


5.67t3-6 


2.07 ±0, 


.63 


66 


0.75 


il5 


0.218 


4.20toi 


25.3ltO 


2.60 ±0.50 


8.80ti;? 


4.94t?-j 


1.90 ±0, 


.51 


66 


0.88 


il6 


0.219 




25.63^^:5 


2.72 ±0.53 


8-99ti;3 


4.53t?:«g 


1.74 ±0, 


.54 


66 


0.72 


il7 
ilS 


0.221 
0.235 


4.10t°J 

4.171°:! 


27.7013-0 
25.03t^;| 


2.83 ±0.32 
2.46 ± 0.34 


10.15+^-^ 
9.74+H 


3.37t4-5 
3.95tJ-f 


1.56 ±0, 
1.83 ±0, 


.33 
.34 


66 
66 


0.92 
0.66 


il9 


0.232 


3.68+^-^ 


30.24tf| 


2.17 ±0.29 


9.0ltJ;° 


4.98t^4 


2.12 ±0. 


.32 


66 


0.91 


i20 


0.259 


3.50t°;i 


30.321^:4 


1.79 ±0.24 


8.42+0-? 




2.02 ±0, 


.26 


66 


1.23 


i21 


0.271 


3-69to3 


27.441^-4 


1.81 ±0.18 


9.86i°| 


4.05tJ;0 


2.02 ±0, 


.20 


66 


0.80 


CM CO 
CM CM 


0.394 
0.379 


4 i4+o'2 
3.68toj 


23.92+2-^ 


2.18 ±0.16 
0.76 ±0.21 


10.94+J-^ 
7.71+H 


2.56t;;? 
3.92+2-^ 


1.22 ±0, 
0.70 ±0, 


.17 

.23 


66 

66 


0.84 
1.22 


i24 


1.940 








6.23t°i 


3.12t°i 


0.19 ±0. 


.03 


68 


0.72 



